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SECTION i 


INTRODUCTION 


An algorithm for the autom-atic synchronization of Omega receivers 
was previously developed and functionally described in Reference 1. 

This algorithm was extended, and its performance for various signal-to> 
noise ratios and station/f requency combinations was determined via 
numerical simulation. Also, the performance of the algorithm was deter- 
mined for receiver-design parameters which effect the algorithm. Param- 
eters considered are the tracking-loop input bandwidth, the data-sampling 
rates, and the length of the data-collecting interval. Values for these 
parameters which offer the best synchronization performancs are suggested 

This report is divided into five main sections. First, the Omega 
navigation system and the requirement for receiver synchronization are 
discussed. Next, the synchronization algorithm is described. Then, the 
numerical simulation and its associated assumptions are described. The 
results of the simulation are presented next. Finally, the suggested 
form of the synchronization algorithm and the suggested receiver-design 
values are presented. Appendix A is a Fortran listing of the synchroniza 
tion algorithm used in the simulation. 
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SECTION 2 


THE OMEOA NAVIGATION SYSTEM 


Omega is a world-wide hyperbolic radio-navigation system. The 
system employs eight stations around the world. The stations are located 
in Norway, Liberia, Hawaii, North Dakota, La Reunion Island, Argentina, 
Trinidad (this is a temporary station and will be replaced with a permanent 
station located somewhere in the South Pacific), and Japan. Each station 
transmits the primary-navigation frequency of 10.2 kHz. The user measures 
the received phase from three or more stations to determine his position. 
The measurement repeats every 360 degrees. Thus, without further informa- 
tion, the determination of position is ambiguous. The 10.2-kHz hyperbolic 
lines-of-position are separated by approximately 8 nmi. Therefore, to 
accurately resolve position ambiguities, the user must have an initial 
knowledge of his position to within >4 nmi. 

The stations also transmit secondary-navigation frequencies of 13.6 
and 11.33 kHz. The secondary-navigation frequencies are used to extend 
the distarce between lines-of-position to 24 and 72 nmi for two- and 
three-frequency receivers, respectively. Thus, to resolve position 
ambiguities, the user must initially only know his position to within 
112 and i36 nmi for two- and three-frequency receivers, respectively. 

To prevent the transmissions from interfering with each other, the 
station/f requency transmissions are time-division multiplexed. Figure 1 
depicts the Omega transmission format. Each station/frequency transmission 
lasts for approximately 1 second. There is a 0.2-second quiet period be- 
tween transmissions. The sequence repeats every 10 seconds. 

Clearly, to naviqate properly, the receiver tracking loops must be 
synchronized to the transmitted format. The transmission format was 
designed such that accurate measurement of the transmission periods of 
any two stations at a single frequency, or of any two frequencies from a 
single station, allow the user to determine unambiguously the starting 
time of the format. When the Omega system is complete, i.e., the South 
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Figure 1. Omega transmission format. 

Pacific station is operational and all the stations are transmitting, 
the starting time of the transmission format can also be unambiguously 
determined with a priori knowledge of the user's approximate position 
and the measurement of a transmission period from a single station on 
a single frequency. Frequently, one or more of the stations is off 
the air for maintenance. Thus, a synchronization algorithm cannot de- 
pend upon a priori knowledge of the stations expected to be received. 
However, the performance of any synchronization algorithm should improve 
with accurate a priori knowledge. 

One method of synchronization is to sense the amplitude of the 
received signal and measure the period of time that it exceeds a j.-re- 
determined threshold. Due to the low signal-to-noise ratios, accurate 
measurements require very narrow bandwidths in the amplitude-senning 
circuitry. In the very low-frequency band, the predominate noise source 
is spike noise ' ' caused by lightning strikes throughout the world. 


* 

Superscript numerals refer to similarly numbered references in the 
List of References. 
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Proper nonlinear diacrimination against spike noise can iir.prove the 
siqnal-tr -noise ratio by 15 dB or more.*^'^^ Thus, many receivers 
employ limiters or clippers in their phase-processing circuitry.**'^* 
Therefore, if the synchronization method requires amplitude measurements, 
both amplitude-sensing and phase-processing circuitry are necessary. 

This increases both receiver complexity and cost. A synchronization 
algorithm that employs phase information does not need the additional 
amplitude-sensing circuitry. Section 3 discusses such an algorithm. 
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SECTION 3 


SYNCHRONIZATION ALGORITHM 


In thl« section, expreeeione for the optimal sine and cosine 
projections of the signal are determined, and two forms of the synchroni- 
zation algorithm are presented. 

3. 1 Op timal Signal-Vector Estimate 

Assume that the signal to be estimated is of the form A sin (wt * 
where A and are constants, and that this signal plus noise is passed 
through a nonlinear system generating an output B sin (u)t * 9) , where B 
is constant and 6 varies with time due to the changing amplitude and 
phase of the noise terms. Figure 2 is a polar representation of two 
such signals. 



Figure 2. 


Polar representation of the Omega s.gnal, A sin 
(u)t i) , and the output of the front end, B sin 
(•-t + 9). 
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Samples B sin (uit *■ 0^) of the output signal are then generated. 
Several such samples are depicted In polar form in Figure ). The 
squared error between each of these samples and th«. Input signal is 

• (A cos ♦ - B cos 6^)^ ♦ (A sin ♦ - B sin 6^)^ (1) 

where is the phase of each sample. Equation (1) yields 

“ (A cos * (B cos - 2AB cos ^ cos 6^ 

*■ (A sin *■ <B sin ~ 2AB sin ^ sin (2) 


The s<]uared error over an Interval in which more than one sample is 
taken is the sum of the individual squared errors, l.e.. 




N N 

NA" cos" t cos^ - 2A cos ♦ B cos 0^ 


,2 __2 


1-1 


i-1 


i-1 


* NA sin 




i-1 


N 

B^ ain^ 0^ - 2A sin ♦ ^ B sin 0^ 

1-1 


(3) 


where N is the number of samples. 



Figure 3. 


Polar representation of the sampled data from the 
front end, B sin (ut 0 ^) . 
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' 1 f fcrcnciatlnq Eq. (3) with respect to A cos | yields 


>c 

jTXcoiTT 


N 

• 2NA cos 3-2 B cos 8^ 
i«l 


( 4 ) 


Settinq this expression equal to zero and solving for A cos 8 gives 
the optimal minimum-squared estimate of A cos i, which is 


N 


A COS 8 


i-i 


B cos 8, 


(5) 


To determine the lea>t-squares estimate of A sin differentiate 
Eq. (3) With respect to A sin t 


3c 

^•(A sin I) 


K 

- 2NA Sin « - 2 ^ B sin 6 
i-1 


i 


( 6 ) 


Thus 


Next, equate this expression with zero and solve for A sin 



A sin $ 


M 




B sin 8^ 


(7) 


3.2 Synchronization Algorithms 

The Omega format consists of eight transmission bursts and eight 
quiet periods which are repeated every 10 seconds. Thus, data from con- 
secutive 10-second segments may be used for synchronization. For example, 
if we call the time at which we start taking data t^ and take 20 seconds 
of data, the samples taken during the intervals (t^, t^ *■ 0.9] and 
It^ 10, t^ * 10.9] may be used to generate the signal -vector estimates 
for Station A. Two methods of employing data from consecutive 10-second 
segments are presented. 

In the first method, data from consecutive 10-second segments is 
used to determine ensemble vector estimates of the received signal at each 
sampling point in the Omega format. Then, a least-square fit of the en- 
semble signal estimates to the local-reference Omega format is performed, 
and the squared error is determined by summing the squared errors between 
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tho ensemble estimates and the local reference, mhe fit is iterated for 
many possible starting times of the local-reference format, and t.he time 
which produces the minimum-squared error is chcn*,; as the correct start- 
ing time. 


The squared error in fitting the signal-vector estimates to the 

local-reference format for the Omega burst, where k • 1 8, 

which correspond to Stations A through H, respectively (see Figure 1), 
is determined by summing the squared errors between the estimates for 
the k^^ burst and the individual sample-point estimates. However, since 
the expression for each of these squared errors will be of a form 
similar to that of Eq. (2), it is not necessary to calculate the ampli- 
tude and phase estimates, only the projections on the sine ind cosine 
axes. These expressions can be determined by expanding Eq. (S) and (7). 
It can be shown that assuming B equals 1, the expanded expressions are 


and 



(A cos 




cos 



(A sin ii), 


N,. r 




w 




sin 6 


j«l L h"l 


hj 


( 8 ) 


(9) 


respectively. Where equals the number of sample points in the k^^ 
burst, W equals the number of consecutive 10-second segments of data 
used, and k is tho same as before. The expressions in brackets on the 
right-hand side of Eq. (8) and (9) are the projections of the j^^ sample- 
point estimatej on the cosine and sine axes, respectively. 

Therefore, from Eq. (2), the squared error for the j*^^ sample point 
in the k^^ Omega burst is 
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- 2 


- 2 




COB 




UOS 


sin 



Thus the squared error 


2 

‘k 


for the k^^ 


Omega burst is 



( 10 ) 



( 11 ) 
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which may b** written ai 



where k, W, j, and h, are the same an before. 

The manner in which the squared error for a quiet period is de- 
fined is not so obvious. Clearly, one may say that during this period 
the signal is by oefinition zero, and therefore the squared error in 
the quiet period, where m • 1, ... 8, is simply proportional to the 
square of the amplitude estimate, i.e. 



where N is the number of sample points (which is the same for all quiet 
periods since they are of equal duration) in the quiet period, and as 
before, W is the number of consecutive 10-second segments employed in 
the computation. 

However, the degree to which zero can be determined in the receiver 
is directly dependent upon the number of samples, and therefore the sampling 
rate. Thus, for a fixed sampling rate there is a Minimum nonzero value 
which will occur during the actual quiet period. Thus, some of the squared 
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error in the quiet periods, as determined by Eq. (13), is due to inherent 
receiver inaccuracies. An alternative method is to assume that the 
squared error in each quiet period is zero, l.e.. 


m 


■ 0 V m 


(14) 


Equations (12) and (13) or (14) may be used to determine the squared 
error for the entire 10-second format, i.e. 

■ Z • E 

k"l m»l 

This process it, iterated for many possible starting times of the 
local-reference format, and the time which generates the minimum-squared 
error is chosen as the correc"; starting time, and the local commutator 
is reset to this starting time . 

In the second method, instead of attempting to determine the 
squared error between the ensemble sample-point estimates and the 
ensemble station estimates (as is done in the first method), the squared 
error between the individual sample points and the ensemble station 
estimates is determined. These errors are summed to produce a total 
squared error for a particular starting time of the local-reference for- 
mat. This process is iterated for many possible starting times of the 
local-reference format, and the time wnich generates the minimum-squared 
error is chosen as the correct starting time. 

Using Eq. (2), (8), and (9), the squared error for the sample 
point in the Omega burst of the h^^ 10-second segment is 



f ^ 

/ w \ ■ 

2 

Sjk ” 


(s Z "hj) 


L j-i 

' h-1 /kJ 


+ (see next page) 
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♦ (co» 


•f 


where 
for the 


’hj’k 


(.in 0^^)^ 


- 2 


- 2 


and W are the same as before. In this case, the squared error 
Omega burst is 


N. 


k W 

k - E L ^hjk 

j-1 h-1 


. W 



/ w 

El 

UL 

j-1 

' h-1 


/ W 

A 

i y 

w 

j-i\ 

> h-1 

! as 

before 

- 

N. W 
k 

•» 

k / 

w 

' 1 i 

w \ W 

Z- 

1 V 

h*l 

'hj’k 

+ N 


cos e 


hj 


).) 


cos e 


hjk 




sin 0 


hjk 


(16) 




j"l \ h«l 


).) 


hj 


'k-< 





^k / w \ 


r ^ w 1 

2 

^ ZU Z 


E E ®hjk 


L j-l\ h-1 'k-i 


L-j-l h-1 J 


r ^k 

/ w 

\ 1 

w 1 


tzi 

L j-1 ' 

J Z 

^ h-1 


E E ®hjk 

.j-1 h-1 

(17) 
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Equation (17) can be rewritten aa 



Equations (18) and (13) or (14) are used to determine the squared error 
for the entire format, i.e. 


■ E 'k • E 

)c*l m«l 

As before, this process is iterated for many possible starting 

times of the local-reference format, and the time which generates the 

minimum squared error is chosen as the correct starting time, and the 

local commutator is reset to this starting time. If Eq. (14) is used 
2 

to estimate e in Eq. (19) , then the minimum Si'uared error will occur 
m 

at the staitinq time which has maximum values for the two right-most 
expressions of Eq. (18). 

The synchronization algorithms discussed previously do not employ 
any type of a priori knowledge. However, the algorithms should have the 
capability of utilizing information concerning the stations expected to 
be received. This information could be manually entered into the receiver 
by the operator. If accurate a priori information is available, the 
performance of the synchronization algorithm should improve. Both of the 
synchronization algorithms discussed can be easily modified to incorporate 
a priori knowledge. Stations whose signals are not expected to be re- 
ceived can be treated in the same manner as a quiet period, and either 
Eq. (13) or (14) can be used to determine the squared errors for those 
stations . 
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SECTION 4 


SIMULATION 


The performance of the eynchronizatlon algorithms presented in 
Section 3 was determined via Monte Carlo techniques. The model used to 
simulate the Inputs to the synchronization algorithms during periods in 
which a signal is received is shown in Figure 4. For each transmission 


R» bM 



N, AND Nj ARE STATISTICALLY INDEPENDENT GAUSSIAN WHITE -NDISE 
PROCESSES WITH ZERO MEANS THEIR VARIANCES ARE OPERATOR SEL- 
ECTABLE A^ AND ARE THE RECEIVED AMPLITUDE AND PHASE FROM 
THE l>"'OMEOA station A^ IS OPERATOR SELECTABLE. AND 0^ IS 
RANDOMLY CHOSEN AT THE START OF EACH HUN 


Figure 4. Model used to simulate the received signal from 
the kth Omega station. 
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burst, the model generates the phase of a constant-amplitude signal from 
a constant-phase constant-amplitude signal and additive white Gaussian 
channel noise. The signal from the burst, sin is broken into 
sine and cosine components, cos sin u.t and sin 6|^ cos u>t, 
respectively. Then, white Gaussian channel-noise processes, cos wt 
and N 2 sin wt, are added to the appropriate quadrature Inputs. The input 
phase to the synchronization algorithms, 6^, is determined by taking the 
arctangent of (Aj^ sin * Nj) divided by (Aj^ cos 6^ * N 2 ) . This produces 
an input phase that is a combination of the phase of tlie received signal 
and a pseudo-random phase due to the noise processes. The amplitude of 
the input signal is set to unity, and hence, the power of the input to 
the synchronization algorithm is fixed, as would be the case if a limiter 
preceded the synchronization algorithm. The amplitude of the received 
signal, A^, is selected by the operator, and may be different for each 
Omega signal. The phase of the received signal, B|^, is randomly chosen 
at the beginning of each Monte Carlo run. is uniformly distributed 

between |-n, it). The input to the synchronization algorithm during the 
quiet periods, and during transmission bursts when no signal is received. 
Is also uniformily distributed between (-it, tt | . 

The input to the synchronization algorithms is sampled every AT 
seconds, where AT inverse is an integral multiple of 10 Hz. For the 
Monte Carlo runs, ATs of 0.1 and 0.05 second were used. The starting 
time of the Omega transmission format can be offset from that of the 
local reference by integral multiples of 0.025 second. 
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SECTION 5 


RESULTS 


This section presents the results of Monte Carlo runs of the simu- 
lation described in Section 4. The performance of the synchronization 
algorithms described in Section 3.2 Is determined for various station/ 
frequency pairs, signal-to-noise ratios, data-sampling rates, and data- 
collection intervals. Two station/frequency combinations were used in 
the comparisons: the A,D 10.2-kllz transmission bursts, and the C,D 

10.2-kHz transmission bursts (the results are valid for any combination 
of station/f requoncies transmitting during these time slots) . Those 
combinations were chosen for tnree reasons. First, they can all be re- 
ceived in the continental United States. Second, the A,D combination 
has a large time separation, 2.7 seconds between transmission bursts, and 
the difference in the duration of the two bursts, 0.3 second, is a maxi- 
mum (see Figure 1). Third, the C,D combination has the minimum separation 
time, 0.2 second, and the minimum difference in burst duration, 0.1 second. 
Thus, synchronization using data from the C,D combination should be- diffi- 
cult, whereas synchronization using data from the A,D combination should be 
relatively easy. Two sampling rates, 0.1 second and 0.05 second, were 
considered. 

Table 1 presents the Monte Carlo results for the synchronization 
algorithm where the squared error is determined from Eq. (12) and (13). 
Notice that the performance of the algorithm increases significantly 
with the use of a priori knowledge. Also, there is a significant per- 
formance difference between the A,D and the C,D combinations when a 
priori knowledge is not used. However, at low signal-to-noise ratios the 
algorithm does not perform very well. 

Table 2 presents the Monte Carlo results for the synchronization 
algorithm employing Eq. (12) and (14). Here again, the performance of 
the algorithm is increased with the use of a priori knowledge, and there 
xs a significant difference in algorithm performance between the A,D and 
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Table 1. Monte Carlo reaulta for the aynchronlzation 
algorithm emnloying Eq. (12) and (13). 


Stations 

AT 

(s) 

No . of ■ 
Periods 

SNR (dB) 

BW - 100 Hz 

No. of Correct Synchronizations 

Without Prior 
Knowledge 
(Out of 50) 

With Prior 
Knowledge 
(Out of 50) 

BB 


10 

0 

50 

50 



10 

-10 

50 

50 



10 

-20 

17 

29 

wgm 


10 

0 

49 

50 


0.1 

10 

-10 

49 

50 

BB 

0.1 

10 

-20 

10 



29 


Table 2. Mont-^. Carlo results for the synchronization 
algorithm employing Eq. (12) and (14). 


Stations 

AT 

(s) 

No. of 
Periods 

SNR (dB) 

BW - 100 Hz 

No. of Correct Synchronizations 

Without Prior 
Knowledge 
(Out of 50) 

With Prior 
Knowledge 
(Out of 50) 

A,D 

0.1 

10 

0 

35 

50 

A,D 

0.1 

10 

-10 

33 

49 

A,D 

0.1 

10 

-20 

6 

7 

C,D 

0.1 

10 

0 

22 

50 

C,D 

0.1 

10 

-10 

19 

50 

C,D 

0.1 

10 

-20 

4 

3 


C,0 combinations without a priori )cnowledge. Comparing Tables 1 and 2 
indicates that the synchronization algorithm used to generate the data 
in Table 1 performs better than the synchronization algorithm used to 
generate the data in Table 2. The only difference between these two 
algorithms is the manner in which they treat quiet periods. The 
algorithm used in the Monte Carlo runs for Table 1 used the data during 
quiet periods to produce an error term; the algorithm associated with 
Table 2 did not. 
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The Monte Carlo result! for the synchronization algorithm employing 
Eq. (18) and (13) are presented in Table 3. In all cases, the performance 
of this algorithm is better, or about the same, as the performance of the 
algorithms associated with Tables 1 and 2. The main difference between 
these algorithms is the manner in which data from consecutive 10-second 
segments is used to generate the squared error. In this algorithm, the 
squared error between the individual sample point and the ensemble 
station estimates is determined. In the algorithms associated with 
Tables 1 and 2, the squared error between the ensemble sample-point esti- 
mates and the ensemble station estimates is determined. 


Table 3. Monte Carlo results for the synchronization 
algorithm employing Eq. (18) and (13). 


Stations 

iT 

(B) 

No. of 
Periods 

SNR (dB) 

BW - 100 Hz 

No. of Correct Synchronizations 

Without Prior 
Knowledge 
(Out of 50) 

With Prior 
Knowledge 
(Out of 50) 

A,D 

0.1 

10 

0 

50 

50 

A,0 

0.1 

10 

-10 

50 

50 

A,D 

0.1 

10 

-20 

18 

35 

C,D 

0.1 

10 

0 

50 

30 

C,D 

0.1 

10 

-10 

49 

50 

C,D 

0.1 

10 

-20 

22 

36 


Table 4 presents the Monte Carlo results for the synchronization 
algorithm employing Eq. (18) and (14). This algorithm differs from 
the algorit)un associated with Table 3 in that it does not make use of 
data taken during quiet periods. The performance of this algorit .n is 
better than that of any of the other algorithms. 

All of the results presented sample the input signal at a 10-Hz 
rate. The effects of a 20-Hz sampling rate for the algorithms associated 
with Table 4 are indicated .in Table 5. From Table 4, the probability of 
synchronization for the A,0 station combination with a sampling rate of 
10 Hz and a signal-to-noisc ratio of -20 dB in a 100-Hz bandwith, is 60%. 
The corresponding probability of synchronization for the C,0 station 
combination is 46%. For a 20-Hz sampling rate and the same signal-to- 
noise ratio, the probability of synchronization for the A,D and C,D station 
combinations was 70% and 60%, respectively. With a priori knowledge, this 
algorithm was able to correctly identify the starting time of the Omega 
transmission format every time. 
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Table 4. Monte Carlo reaults for thv. aynchronlaation 
algorithm employing Eq. (18) and (14). 


Stations 

AT 

(s) 

No. of 
Periods 

‘ SNR (dB) 

BW • 100 Hz 

No. of Correct 

Synchroni rations 

Without Prior 
Knowledge 
(Out of SO) 

With Prior 
Knowledge 
(Out of 50) 

A,D 

0.1 

10 

0 

SO 

50 

i A,D 

0.1 

10 

•10 

>0 

so 

A,0 

0.1 

10 

-20 

30 

46 1 

C,D 

0.1 

10 

0 

50 

50 1 

C,D 

0.1 

10 

-10 

50 

50 1 

C,D 

0.1 

J 

10 

-20 

23 

“ 1 


Tabic 5. Monte Carlo reaulta for the synchronization 
algorithm employing Eq. (18) and (14). 


1 

1 

Stations 

— 

AT 

(s) 

No. of 
Periods 

SNR (dB) 

BW - 100 Hz 

No. of Correct 

Synchronizations 

Without Prior 
Knowledge 
(Out ol 10) 

With Prior 
Knowledge 
(Out of 10) 

A,0 

0.05 

10 

0 

10 

10 

1 A,D 


10 

-10 

10 

10 1 

A,D 

0.05 

10 

-20 

7 

10 

C,D 

0.05 

10 

0 

10 

10 

C,D 

0.05 

10 

-10 

10 

10 

C,D 

0.05 

10 

-20 

6 

10 


As yet, the effects of the time interval over which data is col- 
lected have not been discussed. Table 6 presents the Monte Carlo results 
of the synchronization algorithm associated with Tables 4 and S for three 
data-collection intervals, 100 seconds (which correspond to 10 consecutive 
10-second segments of Omega data), 10 seconds, and 10 seconds. Clearly, 
for low signal-to-noise ratios, 10 or 30 seconds of data is not adequate, 
whereas, with 100 seconds of data, the algorithm performed well. 
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Tabl« 6. Mont* Carlo results for the synchronisation algorithm employ- 
ing Eq. (18) and (14) for various data-collaction intervals. 


Stations 

AT 

(.) 

1 

No. of 
Periods 

SNR (dB) 

BW - 100 Hz 

No. of Correct Synchronizations 

Without Prior 
Knowledge 
(Out of 10) 

With Prior 
Knowledge 
(Out of 10) 

A.D 

o 

o 

10 

-17 

i 

9 i 

10 

A.O 

o 

o 

3 

-17 

2 

7 1 

A.D 

0.05 

1 

-17 

1 : 

5 

A.D 

0.1 

10 

-17 

9 

10 

A.D 

0.1 

3 

-17 

1 j 

1 

6 

A.D 

0.1 

1 

-17 

L_. i. 

0 


Since the algorithm associated with Tables 4 through 6 does not 
ma)(o use of the data coliicted during quiet periods, it is possible for 
a station burst of 1.0-second duration to produce zero error when the 
data is fit to a local-reference burst with a 0. 9-second duration. This 
type ot error can occur when the Omega transmission format is offset 
from the local-reference format by one-half the sampling period. In the 
same manner, station bursts of 1.2 and 1.1 second will produce zero 
error when compared with local-reference bursts at 1.1 and 1.0 second, 
respectively. However, this effect decreases as the sampling period 
decreases, and is less pronounced when the stations received are not con- 
secutive transmission bursts in the Omega format. Table 7 presents the 
Monte Carlo results of this algorithm for a 25-ms offset between the 
transmitted format and the local-reference format. A comparison of 
Tables 5 and 7 indicates that with a sampling period of 0.05 second, the 
effect of the 25-ms offset is negligible. Table 8 presents the Monte Carlo 
results of the synchronization algorithm that employs Eq. (18) and (13), 
and thus uses the data in the quiet periods to generate an error. Clearly, 
the performance of this algsrithm is worse than the performance of the 
algorithm associated with Table 7, which does not make use of data in the 
quiet periods to generate an error. 
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Table 7. Monte Carlo reaulta for the aynchronliatlon algorithm employ- 
ing Eq. (18) and (14) with a 25-ma offaet between the trana- 
mitted format and the local-reference format. 


Stationa 

AT 

(a) 

No. of 
Perioda 

SNR (dB) 

BW • 100 Hz 

No. of Correct 

Synchronize tiona 

Without Prior 
Knowledge 
(Out of 10) 

With Prior 
Knowledge 
(Out of 10) 

] 

A,D 

0.05 

10 

0 

10 

10 

A.D 

0.05 

10 

-10 

10 

10 

A.D 

0.05 

10 

-20 

7 

10 

C.D 

0.05 

10 

0 

10 

10 

C,D 

0.05 

10 

-10 

10 

10 

C.D 

0.05 

^ - 

10 

-20 

5 

10 


Table 8. Monte Carlo reaulta for the aynchronization algorithm employ- 
ing Eq. (18) and (13) with a 25-ma offaet between the trana- 
mitted format and the local-reference format. 


stationa 

AT 

(8) 

No. of 
Periods 

SNR (dB) 

BW - 100 Hz 

No. of Correct 

Without Prior 
Knowledge 
(Out of 10) 

Synchronizations 

wTth Prior 
Knowledge 
(Out of 10) 

A.D 

0.05 

10 

0 

10 

10 

A.D 

0.05 

10 

-10 

10 

10 

A.D 

0.05 

10 

-20 

5 

8 

C.D 

0.05 

10 

0 

10 

10 

C.D 

0.05 

10 

-10 

9 

10 

C.D 

. ^ 

0.05 

i 

10 

-20 

2 

8 
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SECTION 6 


SUGGESTED FORM OF SYNCHRONIZATION ALGORITHM 
AND ASSOCIATED RECEIVER-DESIGN PARAMETERS 


In Section 5, it was shown that for thu eases considered the form 

of the synchronization algorithm with the best performance attempts to 

2 2 

minimize Cq. (19), where f is determined by Eq. (10) and equals zero 
for all m. Equation (18) is minimized when the two right-most expressions 
are maximized. Thus, in terms of performance, an equivalent implementation 
of the algorithm la to choose the starting time of the local reference 
which maximizes the quantity 


8 ( 

t-i ( 


N W 

ZZ 

J«1 h»l 


cos 9 


hj 


(t W 

ZZ 

j-l h-l 


lin 8 


hj 




( 20 ) 


where k, Nj^ , and W are the same as before. This implementation has 
the advantage of requiring less computer operations, and thus less time, 
than the implementation that employs Eq. (18). The results reported in 
Section S indicate that a good choice for W, the number of consecutive 
10-aecond a pnents of data collected, is 10. The results also indicate 
that the input bandwidth to the synchronization algorithm should be 
10 Hz, and the data-sampling period should be 0.05 second. However, if 
a wider bandwidth is desired in the phase-trac)(ing circuitry, then the 
sampling period can be decreased, and consecutive data samples can be 
averaged to generate 0.05-second samples. For example, if a phasc- 
trac)(ing bandwidth of 100 Hz is desired, then the sampling period should 
be 0.005 second, and 10 consecutive data samples should be averaged to 
form a 0.05-second sample. 
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APPENDIX A 

FORTRAN LISTING OF THE SYNCHRONIZATION ALGORITHM 
WITH THE BEST MONTE CARLO RESULTS 
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90002960 
00002970 
00002980 
00002990 
00002900 
00002910 
00002920 
00002930 
00002940 
0000295C 
00002960 
00092970 
00002980 
00002990 
00003000 
00003010 
OOOP 3020 
000030 30 
00003040 
00003050 
00003060 
C0003070 
0CC03090 
O000309O 
C0003100 
00003110 
00003120 
00003139 
00003140 
00093150 
00003160 
00003170 
00003180 
00003190 
00003200 
00C01210 
C0003220 
00003230 
C0003240 
00003250 
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; YAVG«YA»8*Y (K2) 

nsK.-ioArcmi (!2) ) 

XAVr.-YAVS/Rmi 

TAV<5«YA»8/PNi 

AfST (T?) »SOhT(TAVr,**2»YAfG**2) 
? PI':-ST(T2) •ATA*>2(YAVO,XA»(J) 

P7~»»R»I 


?VD 

5ftBhr)«T:tr rcopsv 

r 

« 

R“~tins 

c 

fN 

STBRIOTINE ICOpMD 


C 

RrAL BP(IOPC) 

R7AL ’MCCO) 
p'AL :'a:aip(1''00) 

PL'AL OATPHin?P''» 

CTB'inS /DATBAS/On 
C^niON /EE/S 

CONEON /OATPOL/rA*A*r,PATPHI 
COEEON /IHCCnO/IE 

rU'JT VALENCE (PP(0231) ,ERBPP) , (PB (0 232) , P AT AN) , (DD (02 3 3) , DATPH) 

HHOR-r (IN) 

ATAN«nATAHP (IE) 

\TPH»OATPHI (l.E) 

pr-nRn 


SEP 


C00032B0 
03003270 
00003280 
00003270 
0'' I' 0 3 300 
00003310 
00003320 
00001330 
00003340 

oooonso 

0000 136E 
0C003370 
00003380 
0000 3390 
00003400 
000034TO 
0000 34;’0 
00003430 
P0O0344* 
0000345C 
P000346« 
00003470 
C0003480 
C0003470 
00003500 
00003510 
00001520 
00003530 
00003540 
CC003550 
00003560 
00003570 
00003580 
00003590 
00003600 
00003610 
00003620 
00003630 
0000 3640 
00003650 
00003660 
00003670 
00003680 
00003690 
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